High-energy magnetic excitations in Bi2Sr2CaCu208+(5: Towards a unified description 
of its electronic and magnetic degrees of freedom 
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We investigate the high-energy magnetic excitation spectrum of the high- Tc cuprate superconduc- 
tor Bi2Sr2CaCu208+(5 (Bi-2212) using Cu L3 edge resonant inelastic x-ray scattering (RIXS). Broad, 
dispersive magnetic excitations are observed, with a zone boundary energy of '-SOO meV and a weak 
dependence on doping. These excitations are strikingly similar to the bosons proposed to explain the 
high-energy 'kink' observed in photoemission. A phenomenological calculation of the spin-response, 
based on a parameterization of the ARPES-derived electronic structure and YRZ-quasi-particles, 
provides a reasonable prediction of the energy dispersion of the observed magnetic excitations. These 
results indicate a possible unified framework to reconcile the magnetic and electronic properties of 
the cuprates and we discuss the advantages and disadvantages of such an approach. 



Of all the high- Tc cuprate superconductors, Bi-2212 is 
the most easily cleaved to reveal atomically flat surfaces 
and is the preferred material for many angle-resolved 
photo-emission (ARPES) and scanning tunneling spec- 
troscopy studies of the cuprates. Much is known, there- 
fore, about the electronic structure of this cuprate. In 
particular, ARPES studies have revealed arcs of quasi- 
particle spectral weight at the Fermi energy [IHS], which 
on closer examination have been described as pockets [4] . 
These have strong spectral intensity at the front surface 
(the 'arc') and negligible spectral intensity on the back 
surface [H |5] . This behavior was predicted by the Yang- 
Rice- Zhang (YRZ) ansatz as corresponding to pieces of 
Fermi surface characterized by the (front surface) poles, 
and ('ghost' surface) zeros, of the single particle Green's 
function. |6]. 

In contrast to the electronic structure, however, there 
is scant information about the high-energy magnetic re- 
sponse (> 100 meV) in Bi-2212. This is because Bi- 
2212 is challenging to grow as large single crystals, which 
are required for inelastic neutron scattering experiments. 
Thus only a few studies of the magnetic excitations in 
the cuprates have opted to measure Bi-2212 [SHU] and 
these have focused around the antiferromagnetic order- 
ing wavevector (0.5,0.5), and on energy transfers below 
about 80 meV. The lack of information on the high en- 
ergy magnetic response and its momentum dependence 
away from (0.5,0.5) of Bi-2212 is particularly unfortu- 
nate, given the well characterized electronic structure of 
this cuprate. Such a measurement is also highly rele- 
vant for theoretical proposals for magnetically-mediated 
high- Tc superconductivity. Many such theories require 
the existence of magnetic modes with significant spectral 



weight over a large fraction of the Brillouin zone in the 
optimally doped cuprates [13 . 

The dramatic improvement in recent years of the en- 
ergy resolution of Cu T3 edge resonant inelastic x-ray 
scattering (RIXS) [14-19 means that RIXS now provides 
a useful complement to neutron scattering. It is partic- 
ularly well suited to very small sample sizes < 100 fim 
and can access high energy magnetic excitations above 
100 meV. 

This letter reports the first measurements of the high 
energy magnetic excitations in Bi-2212. We observe 
magnetic excitations with a zone boundary energy of 
~300 meV in both underdoped and optimally doped Bi- 
2212, much like a broadened remnant of the magnon ex- 
citation in undoped cuprates such as La2Cu04 [20) [2T]. 
These excitations are strikingly similar to the boson pro- 
posed to explain the high-energy 'kink' observed in pho- 
toemission [22I [23]. We then take a parameterization 
of the measured electronic structure of Bi-2212 [5] and 
use it to predict the magnetic excitation spectrum in an 
approach based on the Yang- Rice- Zhang ansatz for the 
electronic structure [6 . Such an approach provides a 
reasonable prediction of the energy dispersion of the ob- 
served magnetic excitations. 

We prepared optimally doped Bi-2212 with Tc = 92 K 
(OD92 K) using the floating zone method [25 . These 
crystals were heated under vacuum to produce under- 
doped samples Tc < 2 K (UDO K). Following previous 
work [5], we estimate the hole doping of these samples 
as X = 0.16 and x = 0.03 respectively. Samples pre- 
pared in this way were previously studied using ARPES 
and showed well-defined electronic states, attesting to the 
high sample quality [5 . As a precaution we cleaved the 
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FIG. 1. (Color online) (a) The experimental geometry, (b) A 
typical RIXS spectrum of optimally doped Bi-2212 (OD92 K). 
(c) OD92 K sample at Q = (0.40, 0.00) measured with grazing- 
exit geometry {Oi = 115.6°) and tt or a polarized light, (d) 
Comparison of grazing-exit (Oi = 115.6°) vr-polarized spectrum 
and with grazing-incidence (Oi = 14.4°) cr-polarized spectrum 
for the UDO K sample at Q = (0.40,0.00). We follow pre- 
vious conventions [24] and present spectra normalized to the 
integrated intensity of the dd-excitations. 



samples in-situ in ultra high vacuum ~ 10~^ mbar im- 
mediately before performing the experiments. RIXS ex- 
periments were performed using the SAXES instrument 
at the ADRESS beamline at the Swiss Light Source [26] . 
The x-ray energy was tuned to the peak in the Cu L3 edge 
x-ray fluorescence around 931 eV. The experimental res- 
olution of 150 meV full- width-half- maximum (FWHM), 
and the elastic energy, were determined by scattering 
from disordered graphite. We follow previous conven- 
tions [24] and present spectra normalized such that the 
integrated intensity of the dd-excitations is one. The in- 
plane scattering vector, Q, is denoted using the pseudo- 
tetragonal unit cell a = 6 = 3.82 A and c = 30.8 A, where c 
is normal to the cleaved sample surface. The experimen- 
tal geometry is shown in Fig. [TJa). X-rays are incident 
on the sample surface at Oi and scattered through fixed 
20 = 130°. Oi is varied in order to change Qn - the pro- 
jection of the momentum transfer, Q along [100]. The 
x-ray polarization can be chosen parallel (tt) or perpen- 
dicular (cr) to the horizontal scattering plane. All data 
were taken at T = 15 K. 

Figure [TJb) plots a typical Cu L3 edge RIXS spectrum 
of Bi-2212. We observe a relatively flat continuum of 
excitations up to about 8 eV, which are charge-transfer 



excitations from the Cu 3d to the O 2p states [2T. In 
the 1-3 eV range we see optically forbidden orbital dd- 
excitations where the Cu - y'^ hole is excited into dif- 
ferent orbitals [27]. Below 1 eV, single, and multi-spin- 
flip, excitations are seen. Phonons also contribute below 
~90 meV, along with elastic scattering at eV. 

The peak in the mid-IR portion of the spectrum oc- 
curs around 300 meV, near the Brillouin zone bound- 
ary. This is a comparable energy scale to the mag- 
netic excitations seen in other doped cuprates such as 
La2-:,Sr^Cu04 [l5l[2i[29], and the YBCO and NdBCO 
familes of cuprates [161 ES? suggesting that this peak 
corresponds to single magnon scattering. We test this 
assertion by examining the polarization dependence of 
this peak [3TH33] . In Ref. |32ls formalism the dominant 
term in the intensity of spin-flip scattering can be written 
as / oc \el ' Am,n • where and are the ingoing 
and outgoing x-ray polarizations, respectively, and A^^^ 
is a 3 X 3 tensor . Assuming a Cu^"^ ion in tetragonal 
symmetry, Am,n is zero except for Ai^2 and A2,i. These 
rules work well for describing magnon scattering in in- 
sulating cuprates including La2Cu04 [24 , Sr2Cu02Cl2 
[34 , and NdBaCuOe [16^. In near-grazing-exit geometry 
these rules predict that the magnon should be enhanced 
in TT polarized x-rays and suppressed for a polarized x- 
rays. In Fig. [TJc) the ~ 300 meV peak indeed follows 
this polarization dependence. Further, the elastic line 
at eV is suppressed in the tt channel, consistent with 
expectations for non-resonant charge scattering. For a 
near-grazing incidence geometry the magnon intensity is 
enhanced for a polarized x-rays and Fig. [2|d) shows a 
magnon peak at 300 meV, supporting the proposed con- 
tribution from spin-flip scattering. 

In Fig. [2] we examine the Q-dependence of the OD92 K 
and UDO K spectra along the (0,0) (0.5,0) sym- 
metry direction using a grazing-exit configuration and 
TT polarized x-rays, to maximize the spin-flip scattering 
and minimize the elastic line. The peak is seen to dis- 
perse to high energy near the zone boundary, which fur- 
ther substantiates the assignment of the peak to sin- 
gle spin-flip processes, as multi-spin-flip and phonon 
excitations would be expected to show different dis- 
persions [35l [36]. We find that the peak is signifi- 
cantly broader than the experimental resolution (which is 
150 meV). Thus, although magnon-like excitations per- 
sist in Bi-2212, they are significantly damped. A simi- 
lar paramagnon excitation was also recently reported in 
other doped cuprates Ndi.2Bai.8Cu307, YBa2Cu306.6, 
YBa2Cu408 and YBa2Cu307 M- 

In order to analyze the spectra, we decompose it into 
its different components. The smooth background un- 
derneath the peak, which is largely independent of Q, 
is from charge-transfer processes and the tail of the dd- 
excitations. At zero energy transfer a small amount of 
elastic scattering is visible, although this is strongly sup- 
pressed with TT polarized x-rays. Phonons contribute be- 
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FIG. 2. (Color online) The RIXS spectra of Bi-2212 along the 
(0,0) (0.5,0) symmetry direction (black dots). Left: heav- 
ily underdoped non-superconducting (UDO K); Right: opti- 
mally doped Tc = 92 K (OD92 K). The spectra are decom- 
posed into the dispersive inelastic magnetic scattering (blue 
line), the elastic scattering (red line) and the charge-transfer/ 
(i(i-excitations background (gray dotted line). The data were 
collected at T = 15 K with grazing-exit geometry (Oi > 65°) 
and TT-polarized incident x-rays. 

low ~ 90 meV and multi-magnon processes contribute to 
the high energy tail above the peak [35 . The phonon and 
multi-spin-flip contributions can sometimes be separated 
from the single spin flip contribution in undoped cuprates 
for which the single spin-flip excitation is a well-defined, 
resolution-limited magnon ^5l [161 ttH 121] • In doped 
cuprates, however, the peak broadening makes this diffi- 
cult so, although the peak in scattering can be assigned 
to spin flip scattering, the contributions of phonons and 
multi-spin-flip processes cannot be unambiguously sepa- 
rated. 

We track the peak in the magnetic scattering as de- 
noted by the green arrows in Fig |2] (see the supplemen- 
tary materials for further details) and plot the disper- 
sion in Fig |3] for the (a) UDO K and (b) OD92 K sam- 
ples. The resulting dispersion shows the central result of 
this paper: that high energy spin excitations persist from 
(0,0) (0.4,0) in Bi-2212 and are relatively unchanged 
on increasing the hole concentration from underdoped to 
optimally doped. 

There are a large number of different theoretical ap- 
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FIG. 3. (Color online) Comparison between the peak in the 
RIXS spectra dispersion (solid points) and the theoretical cal- 
culations of S{Q,(jo) based on ARPES-derived parameters [5] 
for (a) underdoped non-superconducting and (b) optimally 
doped Tc = 92 K Bi-2212. The error bars represent in uncer- 
tainty in determining the experimental paramagnon energy, 
(c) The YRZ form for the electronic structure of underdoped 
Bi-2212. 

proaches for calculating the spin response of doped 
cuprates, usually starting with the Hubbard model, the 
t-J model, or models based on stripe correlations. These 
models are then solved by mean field techniques or calcu- 
lations based on small clusters. One difficulty is that few 
of these approaches allow for calculations of the electronic 
and magnetic properties using a single set of parameters. 
For example, it has been argued that local moments con- 
tribute to the spin response at low energies [12 and local 
moment stripe-based theories provide a very natural de- 
scription of the 'hourglass' feature seen in the magnetic 
response near (0.5,0.5) t)ut it is difficult to con- 

nect these theories with the itinerant behavior seen by 
ARPES. 

The YRZ phenomenology [6 , which was inspired by 
calculations on arrays of coupled two leg ladders [38 , 
is advantageous in that it provides a relatively simple 
explanation for ARPES data and an easy way to imple- 
ment parametrization for the electronic structure of the 
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cuprates. It was used in just this way for underdoped 
Bi-2212 |5j. Figure [sjc) plots this electronic structure 
for UDO K Bi-2212 {x = 0.03). It has also proved success- 
ful in describing many other electronic properties of the 
cuprates [7|. 

In Ref. [39 it was argued that calculations based on the 
quasiparticle band structure given by YRZ are also capa- 
ble of describing the magnetic response of the cuprates. 
Here we test this assertion by comparing the magnetic 
response from YRZ, using parameters set by ARPES ex- 
periments, with RIXS data on samples made in the same 
way. The approach taken in Ref. [39] was to set YRZ 
in the framework of earlier calculations by Brinckmann 
and Lee [40 for S'(Q,cj), based on a slave boson mean 
field treatment of the t - J model. The key feature sug- 
gested by the success of YRZ is that one must go beyond 
mean field theory and take into account the binding be- 
tween fermionic 'spinous' and bosonic 'holons'. The ba- 
sic element is a particle-hole 'bubble' diagram of YRZ 
propagators. A geometric series of such diagrams is then 
summed to infinite order leading to an 'RPA' expression 
for the magnetic response. (See the supplementary mate- 
rials for a technical description.) The parameters in this 
phenomenological approach were fixed by fits of the YRZ 
model to ARPES measurements on Bi-2212 crystals [5 . 
The fits assume a pseudogap that decreases linearly with 
doping and vanishes in the overdoped regime, consistent 
with other ARPES studies [41]. These parameters were 
then used to calculate S{Q,uj). We find a broad para- 
magnon at high energies. We also note that a highly 
intense mode arises around Q = (0.5,0.5) with an energy 
of ~ 40 meV [39 , similar to the 'resonance mode' seen 
in neutron scattering [S UTOl [12]. Unfortunately, this is 
above the maximum \Q\ accessible in Cu I/3 edge RIXS. 

Figure [3] compares the results of these calculations to 
the peak in the magnetic response observed for (a) UDO K 
and (b) OD92 K samples. The calculations reproduce the 
energy dispersion, except for a small discrepancy in the 
UDO K sample at high Q. The overall level of agreement 
in the energy of the paramagnon excitation is encourag- 
ing given that the calculations are entirely constrained 
by the electronic structure measurements [5 . 

We note that a comparison of the experimental and 
theoretical widths is difficult to perform in a meaningful 
way, as the width of the measured spectrum is dominated 
by the experimental resolution (150 meV) and the width 
is increased by the contribution from phonons and multi 
spin flip processes. Also, YRZ makes an ansatz for the 
coherent part of the single particle Green's function and 
neglects the incoherent piece. The latter takes the form 
of a large 'hump' in ARPES data suggesting that in fact 
lifetime effects play an important role. A more sophisti- 
cated theory, in which this effect is accounted for, would 
probably predict larger widths. 

Experimentally, we find similar intensities in the 
UDO K and OD92 K samples, although it is difficult to 



measure RIXS on an absolute intensity scale. Neutron 
scattering measurements observe a strong suppression 
of spectral weight around (0.5,0.5) in doped cuprates 
[42] . Taken together the x-ray and neutron results then 
suggest a strong Q-dependence to the doping depen- 
dence, with spectral weight largely conserved, except 
around (0.5,0.5). Such a Q-dependence can be repro- 
duced within YRZ ^ , and other theoretical approaches 
such as Ref. [43] . Upon increasing the doping from heav- 
ily underdoped to optimally doped, the calculations pre- 
dict a small drop of 15% in the intensity of the param- 
agnon. In future theoretical work we intend to include 
the possible appearance of an ordered moment in the 
YRZ description at low doping. This should improve the 
accuracy of the theory for Q (0,0), which is currently 
dominated by the large pseudogap. 

Our results also have relevance to ARPES studies. Sev- 
eral measurements of Bi-2212, and other cuprates, ob- 
serve a kink of the electronic bands in the mid-IR en- 
ergy scale [H [231 [44]. In Bi-2212 this feature occurs 
around 340 meV and it is weakly doping dependent [22] . 
These characteristics are similar to the magnetic exci- 
tations near the Brillouin zone boundary reported here, 
and indeed, calculations [43] have shown that a param- 
agnon with just the characteristics of the mode we mea- 
sure would produce the observed phenomenology of the 
high energy kink. However, it remains unclear whether 
the high energy kink is due to electron-boson coupling, 
or to other effects [45] . 

In conclusion, we have used RIXS to determine the 
high energy magnetic excitation spectrum of Bi-2212. We 
find broad dispersive paramagnons with a zone boundary 
energy of ~300 meV and a weak dependence on doping. 
The observed energy dispersion of these excitations are 
reasonably well reproduced by calculations based on the 
YRZ model [6 with parameters fixed by the measured 
electronic structure of Bi-2212 [5 . This indicates that 
this is a promising path towards a possible unified de- 
scription of the magnetic and electronic properties of the 
cuprates, as opposed to models that focus on just one of 
these degrees of freedom and we discuss the advantages 
and disadvantages of such an approach. Future work 
should improve the limitations of the current theory and 
test whether other spectroscopies can be incorporated 
into the same unified picture. 
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Here we provide an additional diagram depicting the 
scattering geometry, a more detailed description of the 
fitting of the RIXS spectra and further details of our pa- 
rameterization of the Bi2Sr2CaCu208+5 electronic struc- 
ture according to the Yang- Rice- Zhang ansatz (YRZ) [1]. 
Then we go on to describe how this provides a basis for 
calculating the magnetic response. 

EXPERIMENTAL GEOMETRY 

Figure 1 shows further details of the scattering geom- 
etry used in the experiment. In addition to the param- 
eters defined in the main text, here we define the total 
momentum transfer Q, which is resolved into (5||, parallel 
to a'^, and Q^, parallel to c'^, J is the angle between Q 
and Q^. In this diagram 7r-incident x-rays have polar- 
ization in the plane of the diagram perpendicular to ki 
and cr-incident x-rays have polarization perpendicular to 
the plane of the diagram. The ADRESS beamline used 
for these measurements employs an APPLE-II type he- 
lical undulator which produces highly (~ 99%) polarized 
x-rays. 

FITTING THE RIXS SPECTRA 

As described in several studies of the cuprates (see for 
example Refs. [2-6]) several different types of scattering 
processes are present in Cu L3 edge RIXS spectra. These 
include elastic processes from specular and diffuse scat- 
tering, phonons below ~ 90 meV, single spin-flip scat- 
tering, multi spin-flip scattering, charge scattering and 
orbital scattering. All the peaks used in the fitting here 
and in previous work [4] were convolved with a Gaussian 
function, which is an excellent approximation to the res- 
olution function of the SAXES instrument used in this 
study. By definition, the elastic scattering is a ^-function 
at eV energy transfer, which leads to a Gaussian when 
convolved with the resolution function. The electronic 
and orbital scattering causes a tail of scattering in the 
low energy region, which can be modeled by a smooth 
function, and in this case we used a third order poly- 
nomial. In YEGG (Ref. 4) the RIXS spectrum could 
be fit with a resolution-limited Gaussian to account for 
the elastic line, an anti-symmetrized Lorentzian to ac- 
count for the phonon, single spin-flip paramagnon and 
multi spin-flip scattering and a smooth background. In 
Bi-2212 we find that the low energy inelastic scatter- 
ing cannot be adequately accounted for with an anti- 
symmetrized Lorentzian as the high energy tail of the 



Vectors in the scattering plane 

A 




FIG. 1. (Color online) The scattering geometry in the ac- 
plane showing the definitions of the different angles and vec- 
tors used in the experiment, ki and kf are the initial and 
final scattering vectors of the x-rays. 20 = 130° is the angle 
between ki and kf. Oi is the angle between the sample surface 
and ki. a* and c* are the reciprocal lattice vectors. Q is the 
total momentum transfer, which is resolved into Q\\, parallel 
to a* , and Qi, parallel to c'^. 



paramagnon extends out too far in energy. Similarly, 
a measurement of Lai.92Sro.o8Cu04 also could not ac- 
count for the low energy inelastic scattering in terms of 
a single anti-symmetrized Lorentzian [3]. Given that the 
phonon, single spin-flip and multi spin-flip scattering in 
Bi-2212 could not be unambiguously separated, we used 
an anti-symmetrized Lorentzian plus two Gaussians to 
provide a smooth line, which accounts for the MIR re- 
gion of the RIXS spectrum. In essence, however, this 
is a phenomenological method to provide a smooth line 
within the confidence limits set by twice the errorbars. 
The paramagnon energy is taken to be the peak of this 
smooth line. The error is estimated by adding the un- 
certainty in the fit procedure to the uncertainty in de- 
termining the zero energy transfer energy via scattering 
from carbon tape. 
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YRZ PARAMETERIZATION OF THE 
ELECTRONIC STRUCTURE OF Bi-2212 

In the YRZ model the coherent single particle excita- 
tions are described by a Green's function for spinons 



TABLE L The bare parameters. 



to 



J 



Ao 



3J/2.5 -0.3to 0.2to 0.12 eV 0.338 0.3to 



G{uj, k) 



(1) 



where ^(/g) is a hopping dispersion (including terms up 
to next- next- nearest neighbor) renormalized by doping 

m=Uk)+^'{k), (2a) 
^o(^) = -2t{x){coskx cos ky), (2b) 
= -4:t\x) cos kx cos ky 

-2t" {x){cos2kx + cos2ky) - fip^ (2c) 

and 

I]fi(a;,fc) = |Afi(fc)|2/[a; + Co(fc)], (3) 

is a self energy term that depends on the pseudogap 

Ai^(/g) = Ao(x)(cosA;^ -cosky). (4) 

The doping dependence is included via Gutzwiller factors 
[7], 



9t 



2x 

1 + x 

4 

(i + xy 



(5a) 
(5b) 



which renormalize the bare hopping and exchange pa- 
rameters to 5 5 ^0 



t{x) =gtto + -gsJx^ 
t'{x) =gtto, 



(5c) 

(5d) 
(5e) 



with X = 0.338 [8]. The bare parameters arise from an 
underlying t - J Hamiltonian description 

Ht-j = Ps{-Z to,ijcUjs + J Z 5i • Sj)Ps, (6) 

ijs (ij) 

where Pg is a projection operator onto the singly occupied 
subspace. Note that the exchange J enters the dispersion 
as a hopping-like term due to a mean field factorization 
of the spin interaction [1] . 

The precise form of the pseudogap in the cuprates is 
still an issue of debate, especially the question of whether 
the pseudogap closes at optimal doping or in the over- 
doped regime where superconductivity disappears [9, 10]. 
The original formulation of YRZ assumed the former [1] 
but ARPES measurements tend to side with the latter 
scenario [10] and indeed this is the finding of ARPES 
measurements on Bi-2212 samples made in the same way 



as the sample we study here [11]. Given that our aim 
is to reconcile ARPES measurements of the electronic 
structure of Bi-2212 with our RIXS measurements of the 
magnetic excitations, we follow the convention of Ref. 
[11]. Namely we choose a simple linearly decreasing phe- 
nomenological form for the pseudogap 



Ao(x) = Ao(l -x/xcrit), 



(7) 



which disappears at Xcrit = 0.20, the overdoped edge of 
the superconducting dome [12]. The overall magnitude, 
Ao, is chosen by fitting to the measured electronic struc- 
ture [11] as explained later. We note that our results are 
not highly sensitive to the precise value of Xcrit- 

ARPES, when corrected for matrix element effects and 
divided by the Fermi occupation factor, measures the 
spectral function A{k^uj)^ which is directly related to the 
Green's function as A{k^uj) = -(l/7r)ImG(/g, cj). This re- 
lationship was used in [11] to determine the parameters 
for the YRZ model applied to Bi2Sr2CaCu208+5. Table I 
shows this bare parameter set [12]. As emphasized in the 
text, we use the doping values obtained in the fitting of 
the ARPES data: For the underdoped sample x = 0.03 
and at optimal doping x = 0.16. The chemical potential, 
fip is set according to Luttinger's sum rule [1] self con- 
sistently and found to be -0.036 and -0.054 for x = 0.03 
and 0.16 respectively. We plot the resulting YRZ forms 
for the electronic structure of Bi-2212 in Fig. 2. 

Unfortunately, the YRZ ansatz in its original form is a 
single particle Green's function and does not constitute 
a starting point for calculating the magnetic response. 
However, a slave boson treatment of the t - J model, in 
which fermionic 'spinous' are bound to bosonic 'holons' 
reproduces the YRZ propagator [13]. This provides a 
connection to standard many body theoretical techniques 
[14, 15]: one can calculate S{Q,lj) via a resummation of 
particle-hole bubble diagrams [14]. 

The YRZ ansatz factors into two effective bands and 
can be extended to include superconductivity by treat- 
ing the system as a two band superconductor [16]. We 
write the superconducting gap in the i band as As^i{k) = 
As^i(x)[cos(kx)-cos(ky)]. Because the upper, (+), YRZ 
band is far above the Fermi energy we set As^+ = and 
take As - = 0.07to(l - 91{x - Xcrit)^) in the supercon- 
ducting state [11, 12]. In the heavily underdoped, non- 
superconducting case As^+ = Ag- = 0. The propagator 
is then written as a Nambu-Gor'kov (matrix) Green's 
function (with momentum labels suppressed and Pauli 



FIG. 2. (Color online) The YRZ form for the electronic structure of Bi-2212, which produces two effective bands with the 
Fermi energy lying in the lower band. The resulting Fermi surface consists of hole pockets, with an area proportional to doping. 
The intensity of color is proportional to the coherent quasi-particle weight. The lower, blue band has strong quasi-particle 
weight on the side near (0, 0) and weak quasi-particle weight on the far side towards (0.5, 0.5). The gray surface represents the 
Fermi energy, which forms the Fermi surface outlined in yellow. The gap between the upper and lower bands closes along the 
Brillouin zone diagonal (0,0) (0.5,0.5) at the top of the pocket. A key feature of the YRZ model is that the quasiparticle 
weight vanishes at the back of the pockets. The left hand plot shows the bands for the underdoped, x = 0.03 case and the right 
hand plot shows the optimally, x = 0.16, doped case. 



matrices r) 



with 



u;, = ^±^^"(fe)2 + A|j(fe), 
2 2^akr + Alik) 



(9) 



and we have defined ^{k) = [^{k) + ^o{k)]/2. 

The lowest order contribution to the magnetic response 
at T = is 

dk ^ Zi(k)zj(k+ Q) 



Ji(k)ujj(k+ Q) + As,i(k)Asj(k+ Q) - 
Es,i(k)Es,j(k+ Q) 

<[ I 

\-uj-E,j{k+ Q)-E,,i{k) 



-i ] 

+ Q) + E.,j(k)y 



< + E,j{k+Q) + E,,i{k) 



(10) 



where 



E,,{k) =^coUk) + Al^{k). 



(11) 



On setting = one recovers the T = normal state 
particle-hole bubble diagram, composed of YRZ propa- 
gators. We evaluate Ils{^^ Q) numerically by summing 
over a grid of 1000 x 1000 points in k space. 

The dynamical structure factor is given by the RPA 
resummation 



7r»?-o 1- J(Q)n*'(w + i7?, Q)' 



(12) 



with an antiferromagnetic exchange J( Q) = 2 J(cos Qx + 
cos Qy). Here we write J to indicate that the bare ex- 
change J that appears in Eq. 6 should be renormalized 
by higher order terms (such as interactions between holes 
and spin excitations that disrupt antiferromagnetic cor- 
relations) [14]. The neglect of these higher order terms 
is a limitation of our approach. In fact in Ref. [14] the 
necessary renormalization factor was found to be quite 
large, ~ 1/3, in order to produce the transition to antifer- 
romagnetic order at a low doping. Treating J as a fitting 



4 



parameter for the ~ 40 meV resonance measured at op- 
timal doping {x = 0.16) we find J = 135 meV, compared 
to J = 120 meV. Further discussion of these issues can 
be found in Ref. [13]. For the energies and wave vectors 
relevant to the RIXS data presented in this work the re- 
summation is negligible. It does however have some effect 
at low energies near (0,0), causing a small broadening 
in energy. Inserting the numerically evaluated bubble, 
Eq. 10, into Eq. 12 we arrive at the theoretical results for 
the magnetic response shown in Fig 3. of the letter. 
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